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. Abstract. I review some recent cosmological studies based on redsliift 

' surveys of radio sources selected at low frequencies. The accretion rate 
onto the central black hole is identified as the basis of a crude physi- 

^ I cal division of the low-frequency population into two sub-populations: 

^ ■ the first comprises 'Eddington-tuned' (high accretion rate) quasars and 

I their torus-hidden counterparts; the second comprises 'starved quasars' 

. like M87. There exist remarkable similarities between the shapes and 

I I evolutionary behaviours of the luminosity functions of radio sources and 

^ ■ radio-quiet quasars; all luminous AGN are ~ 300— times rarer at epochs 

r — I corresponding to z ~ than at z ~ 2.5. I argue that any evidence that 

■ quasars were intrinsically rarer at z ~ 5 than at z ~ 2.5 is as yet both 
I tentative and indirect. A simple calculation suggests that the radio source 

Q ■ population has been over-looked as a potentially important contributor 

I to the entropy budget of the Universe. A recent sub- mm survey of radio 

■ sources is used to demonstrate a connection between the events which 
I trigger jets and intense bursts of star formation, and a close link between 

CLc the histories of star formation and AGN activity is proposed. I discuss 

Q ! the aims and methods of future large redshift surveys of radio sources, 

i-{ ' emphasising the importance of dovetailing these with the development 

c/3 . of robust physical models for radio sources and their epoch-dependent 

^ ' environments. 
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I 1. Redshift surveys of low- frequency radio sources 

We were all excited at this meeting by the prospect of the GMRT pushing low- 
frequency astronomy to sensitivities significantly below those of the previous 
generation of surveys: e.g. the ~ 5a detection levels of 151 MHz flux density 
5151 ~ 100 mJy reached by the Cambridge 7C survey (e.g. McGilchrist et al. 
1990). Nevertheless, the next generation of low-frequency continuum surveys will 
continue to be dominated by extragalactic sources whose large radio luminosities 
come ultimately from energy-releasing processes associated with super-massive 
black holes. In this sense most of the sources in a GMRT survey will not be so dif- 
ferent from those in the 3CRR survey and, as we have recently been reminded by 
Longair (1999), these sources have been studied for nearly forty years. So, what 
new things can we learn from low-frequency radio surveys deeper than 3CRR? 
As a partial answer to this question I will develop two themes illustrated by 
the increased coverage of the 151 MHz luminosity^Lisi, redshift z plane (Fig. 1) 
provided by recently completed redshift surveys of radio sources. 
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Figure 1. The coverage of the 151 MHz luminosity L151 versus red- 
shift z plane for existing, planned and possible future redshift surveys 
of low-frequency selected radio sources. The number of objects in each 
survey is given in brackets, slightly above the locus of the relevant 
5151 lower limit (taking a = 0.8). References for the named surveys 
are: 3CRR (triangles) - Laing, Riley & Longair (1983); 6CE (circles) 
- Rawhngs, Bales k Lacy (2000); 7CRS (shaded region) - Willott et 
al. (2000d), Lacy et al. (1999), Blundell et al. (2000); TOOT - Sec. 6. 

The first theme concerns the demography of the low-frequency radio pop- 
ulation. The tight distribution of the 3CRR points on Fig. 1 results from the 
conspiracy of three features of the radio luminosity function (RLF) — the den- 
sity PR = [d/9(Li5i, 2:)/dlog^o ^151] of radio sources per unit co-moving volume 
and per unit logarithm of radio luminosity — namely: (i) the points hug the lo- 
cus of the 3CRR 5151 limit because at any given z, the RLF is a steeply declining 
function of L151 ; and (ii) the points are fairly uniformally distributed along this 
locus because out to z ~ 2 both the normalisation of the RLF and the observable 
co-moving volume are steeply increasing functions of z. A corollary of these fea- 
tures is that the RLF is well determined close to the locus of the Si^i limit but 
estimates elsewhere in the Li^i-z plane require extrapolations of what is clearly 



^Throughout this paper I take Ho = 50 km s~^Mpc~^, 0,m = 0.3 and Qa = 0.7; radio spectral 
index, a, is defined by Si, oc where Si, is the flux-density at frequency v; and radio 

and optical narrow- line/bolometric luminosities are measured in units of W Hz~^ sr"'^ and W 
respectively. 
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a rapidly varying function (although data like the low-frequency radio source 
counts provide additional constraints) . Constraining the RLF at regions of the 
-^151-2 plane corresponding to higher aerial densities of radio sources (e.g. the 
7C Redshift Survey, 7CRS, region of Fig. 1) is essential for many studies which 
utilise radio sources as cosmological probes: e.g. gravitational lensing studies 
which constrain cosmic geometry (e.g. Kochanek 1996) and clustering studies 
which map out large-scale structure (e.g. Magliocchetti et al. 1998). 

Of more direct concern to us here is the ability to understand the RLF in 
terms of physical models for radio sources and their cosmic evolution. Ignoring 
for now any variations in radio source environment, what we would really like to 
measure is the density pp{Q, t, z), where Q is the jet power (Scheuer 1974; Rawl- 
ings &; Saunders 1991), and t is the source age, the time after the jet-triggering 
event. Prom pp (Q, t, z) it is possible to calculate physically meaningful quantities 
like the birth function - the number of jets triggered per co-moving volume per 
unit cosmic time — and the total entropy injected into the intracluster medium 
(ICM) by the radio source population. Successful mapping between the 'phys- 
ical' density pp{Q,t,z) and the 'observable' density po^Li^i, D, z) (and hence, 
by marginalising the projected radio source size D, pr) requires both a good 
physical model and samples providing adequate coverage of Fig. 1. Inadequate 
coverage leads to unbreakable degeneracies between the key parameters: Blun- 
dell &: Rawlings (1999) have recently emphasised one such effect which means 
that the known highest-z radio galaxies are biased towards objects with ex- 
tremely high Q and relatively small t, some complicating effects of which will 
be discussed in Sec. 5. 

This introduces the second important theme illustrated by Fig. 1 - the use 
of radio sources as probes of the very high redshift Universe. The highest redshift 
object in the 3CRR sample is at 5; 2. The lack of higher-z 3CRR objects is not 
due to an abrupt high- 2; cut-off in the RLF, but due instead to the intrinsic rarity 
of radio sources with the extreme radio luminosities [log]^o(-^l5i)~29] required to 
exceed the S151 limit of 3CRR at high redshift. At a given redshift the finite 
(co-moving) volume on our light cone depends only on the assumed cosmological 
model, and since the 3CRR survey covers a good fraction of the entire sky, the 
only way of finding significant numbers of higher-z radio sources is by studying 
objects with luminosities similar to those of the 3CRR objects. This requires 
spectroscopic surveys of fainter radio-selected samples; the Oxford-led 6CE and 
7CRS redshift surveys (see Fig. 1) have demonstrated that one can achieve 
close to 100% redshift completenesses for such surveys using 4-m class optical 
telescopes. This has allowed the construction of sizeable and complete samples 
of 2; > 2 radio sources from follow-up of fairly small sky areas. 

2. Sub-populations of low-frequency radio sources 

The traditional way of dividing the low-frequency radio population into two 
distinct sub-populations is via the Fanaroff &: Riley (1974) scheme based on 
radio structure: crudely, FRII galaxies have jet-terminating compact hotspots 
at their outer extremities, whilst FRI galaxies do not. Fanaroff &: Riley showed 
that this FRI/FRII structural division corresponds closely to a division in ra- 
dio luminosity at logio(^i5i) ~ ^5.5 although more recent studies (Ledlow & 
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Owen 1996) have suggested an additional dependence on host galaxy luminos- 
ity. However, following Hine & Longair (1979) and Laing et al. (1994), there is 
another, arguably more fundamental, way of dividing the population into two 
classes. I will adopt two underlying assumptions: (i) that all radio galaxies 
harbour super-massive 10^ M©) black holes with a low spread in black hole 
mass (see Willott et al. 1999), and (ii) that higher accretion rates onto the black 
hole boosts: (a) the quasar optical luminosity, and hence the narrow emission 
line luminosity; and (b) the jet power, and hence the radio hnninosity (Rawlings 
& Saunders 1991; Serjeant et al. 1998). I will neglect the effects of scatter in 
correlations between L151 and, for example, the [Oil] emission line luminosity 
L[Oii] (e-g- Baum & Heckman 1989; Willott et al. 1999). 

Above logio(-i^[o//]) ~ 35 [or logio(-i'i5i) ~ 26.5], the 'quasar fraction' - 
the fraction of objects with observed broad emission lines - is 0.4 regard- 
less of luminosity or redshift, a result which is consistent with a simple unified 
scheme in which each object has an obscuring torus with a half-opening angle 
^trans ~ 53° (Willott et al. 2000b). In other words this appears to be a homoge- 
neous 'Eddington-tuned' sub-population comprising objects whose black holes 
are accreting at some sizeable and roughly fixed fraction of their Eddington- 
limited rates, and whose nuclei have well-developed obscuring tori. 

Below logiQ{L[on]) ~ 35 [or log]^o(^i5i) ~ 26.5], the quasar fraction drops 
abruptly, and a second sub-population dominates. The archetype of this sub- 
population is the nearby FRI radio galaxy M87 which seems certain to contain a 
super-massive black hole accreting at a tiny (< 10^^) fraction of its Eddington- 
limited rate (e.g. Ford et al. 1994). Chiaberge, Capetti &: Celotti (1999) discuss 
how the seemingly clear view of optical synchrotron from the base of the M87 
jet argues against the existence of an obscuring torus in this radio galaxy. 

Because the break between these 'starved quasar' sources and the Eddington- 
tuned population lies about one dex higher in L151 than the FRI/FRII division, 
there are starved quasars with both FRI and FRII radio stmcturc, whereas the 
Eddington-tuned sources are almost exclusively FRIIs. The existence of FRIIs 
lacking an intrinsically bright quasar nucleus is supported by cases in which nar- 
row emission lines are either weak or absent (Laing et al. 1994), and, as for M87, 
the detection of optical synchrotron from some FRII nuclei rules out some types 
of obscuring tori in some FRII objects (Chiaberge, Capetti & Celotti 2000). 

Being based on accretion rate, the division proposed here is much more 
intimately tied to physical processes near the black hole than the FRI/FRII 
divide which is likely to be a larger-scale magnetohydrodynamical effect driven 
by the dependence of the stability of jets on jet power and environment (e.g. 
Kaiser & Alexander 1997). Since Eddington-tuned radio sources always contain 
powerful, and hence stable, jets it is no surprise that they almost exclusively 
develop FRII radio structures. The less powerful jets in 'starved quasars' can 
produce cither FRI or FRII radio sources depending on their precise jet power 
and environment. 

The accretion rate based division proposed here is undoubtedly too simplis- 
tic: as suggested by Rawlings &: Saunders (1991), variations in accretion rate 
from object-to-object might be a far smoother affair than a clean dichotomy. 
Indeed, the possibility that the low quasar fraction at low luminosities is due to 
lower values of ^trans for weaker quasar nuclei, rather than the intrinsic lack of 
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such nuclei, has received support from the recent evidence for accretion-disc ex- 
cited broad hnes in the pole-on low-luminosity radio source BL Lacertae (Corbett 
et al. 2000). However, as as we shall see in Sec. 3, the two sub-population hy- 
pothesis proposed here provides a simple luminosity-dependent density evolution 
(LDDE) model on which to base our first stabs at mapping and understanding 
the cosmic evolution of the radio source population. 



3. Cosmic evolution between z ^ and z ~ 2.5 




Figure 2. Comparison of the himinosity functions of active galaxies 
selected at X-ray and radio wavelengths at z = (left), and z = 2.5 
(right). The upper solid line is a simple scaling of the LDDE2 model of 
Miyaji et al. (2000) for the soft X-ray luminosity function (upper dotted 
line) , adopting a scaling factor ~ 4 to correct for the large population of 
Compton-thin, photoelectrically-absorbed quasars required to provide 
the X-ray background (e.g. Wilman & Fabian 1999); I have assumed 
that « 2 per cent of the bolometric luminosity Lboi of a naked quasar 
emerges in the 0.5 — 2 keV band (Elvis et al. 1994). The lower solid 
line shows the model stccp-spcctrum RLF from Willott et al. (2000c), 
and the lower dashed line is a model RLF for just the steep-spectrum 
quasars from Willott et al. (1998); L151 and dp/dlog^Q L151 have been 
converted to Lboi and dp/dlog^g ^boi respectively, assuming that Lboi = 
1016-2 X L^sf (Willott et al. 1999). 



Willott et al. (2000c) have used the 3CRR, 6CE and 7CRS datasets il- 
lustrated in Fig. 1 (alongside other pertinent data like the low-frequency radio 
source counts) to constrain parametric models for the cosmic evolution of the 
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RLF; these models were built on the assumption that the two sub-populations 
described in Sec. 2 undergo differential density evolution. They found that the 
integrated co-moving number density T = J[dp/dlogiQ Li^i] dlog^gLisi of the 
Eddington-tuned sources rises by about 2.5 dex between z ^ and z ~ 2.5, 
and that the T of the starved quasar sub-population rises by about 1 dex be- 
tween 2; ~ and z ~ 1. The superposition of these sub-populations produces 
an overall steep-spectrum RLF in fair agreement with a previous determination 
based on high-frequency (2.7 GHz) selected samples (Dunlop &: Peacock 1990). 
The 'break luminosity' in the overall RLF is the swap-over point between the 
two sub-populations. Because these sub-populations undergo differential density 
evolution, the break point moves with redshift in a way which is similar to the 
'pure luminosity evolution' (PLE) model fitted by Dunlop & Peacock to their 
data. The Willott et al. LDDE model is more attractive than a PLE model 
because it is hard to see how the latter can provide a physical description of 
evolution given that the ages of FRII radio sources are known (e.g. Blundell &: 
Rawlings 2000) to be a tiny fraction of the Hubble time. 

I have chosen to illustrate these results in a way (Fig. 2) which emphasises 
some remarkable similarities between the evolution of the radio population, and 
quasars in general. First, the shapes of the radio and X-ray luminosity functions 
are almost identical, and it is interesting in this regard that the most recent 
estimates of the X-ray luminosity function (Miyaji, Hasinger & Schmidt 2000) 
prefer LDDE to PLE models. Second, the evolutions with redshift are very 
similar. We may soon become persuaded of a third similarity related to the 
luminosity-dependent obscuration of the quasar population. Studies of the RLF 
(Willott et al. 2000b) show that the fraction of objects with observed broad lines 
rises from ^0.1 at log;^o(^6oi) 39 (the Eddington-limit of a 10^ Mq black 
hole) to « 0.4 at high luminosities. For the X-ray population at logio(-^>6oi) 39, 
the fraction of naked quasars can only be <0.2 since the hard X-ray background is 
dominated by a photoclcctrically-absorbcd, but Compton-thin population (e.g. 
Wilman &; Fabian 1999). Follow-up of Chandra surveys should confirm whether 
or not the fraction of naked quasars rises with luminosity, as is the case for the 
radio population, and as is predicted by the popular 'receding-torus' model 
(Lawrence 1991). First results from Chandra (e.g. Mushotzky et al. 2000; 
Crawford et al. 2000) suggest that it may well prove just as difficult to measure 
redshifts for the X-ray selected 'quasar lis' as has proved the case for their radio- 
loud counterparts in the 7CRS (Willott, Rawlings & Blundell 2000a): nuclear 
obscuration means that quasar lis can masquerade as extremely-red galaxies 
with only weak observable narrow emission lines in the near-IR, especially in 
the key 1.2 < z < 1.7 redshift range. 

One open question concerns the size of any Compton-thick population of 
AGN. High-energy X-ray studies of local AGN (e.g. Maiolino et al. 1995) suggest 
that, at least at low values of L^oi, most nuclei are obscured by Compton-thick 
material, so the normalisation of the quasar luminosity function in Fig. 2 is still 
strictly only a lower limit. Chandra follow-up of AGN selected at low radio 
frequencies might prove vital here since the large radio emitting lobes should 
be unobscured, allowing direct measurement of the distribution in obscuring 
column densities for a randomly-oriented population of quasars. 
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Fig. 2 also illustrates the key remaining challenge for understanding the 
quasar population: why are there seemingly ^ 100— times more radio-quiet 
quasars than radio-loud quasars^? The most naive explanation - that 'radio- 
quiet' phases of quasar activity last ~ 100— times longer than 'radio-loud' phases 

- is not really tenable. Since jet-fed powerful radio sources persist alongside nu- 
clear optical/X-ray quasar activity for timescales of ~ 10^ years (e.g. Blundell, 
Rawlings & Willott 1999; Blundell & Rawlings 2000), then any such radio-quiet 
phase would need to last ~ 10 Gyr, leading to just one generation of quasars 
and a well-known problem explaining the mass density of super-massive black 
holes in the local Universe (e.g. Richstone et al. 1998). The cosmic evolutions of 
the radio-loud and radio-quiet populations (Fig. 2) seem to be driven by a com- 
mon cosmic cause, presumably ~ 10^-yr timescale variations in the trigger rate 
convolved with shorter-timescale (~ 10^ yr) effects due to finite quasar lifetimes. 

A more likely explanation for the large fraction of radio-quiet to radio- 
loud objects follows from the invocation of a minimum black hole mass for the 
development of powerful jets. Consider a toy model in which quasars divide 
cleanly between radio-quiet quasars and radio-loud quasars at a critical black 
hole mass Mbh ~ 10^'^ M©: because of the steepness of the mass function for 
black holes (e.g. Salucci et al. 1999), for every radio-loud quasar there would be 
a factor >100-times more radio-quiet quasars with black hole masses in the range 
8.5 < log^Q Mbh < 9.5. Thus, a spread of only ~ 1 dex in the ratio of Lboi/MeH 
would cause the radio-quiet quasars to swamp radio-loud quasars even at the 
highest values of Ly,oi (see also McLure et al. 1999). There is some indirect 
support for the central idea here - that, compared to radio-quiet quasars, the 
distribution in Mbh for radio-loud quasars has a high mean and a low variance 

- from the normalisation and spread in the near-infrared Hubble diagram for 
radio galaxies (e.g. Bales et al. 1997). Although it seems hard to envisage a 
critical physical switch allowing powerful jet activity only at Mbh~10^'^ Mq, 
it may appear to be present because of the special formation history of the 
most massive black holes. As explored by Wilson & Corbett (1995), the major 
merger of two massive dark haloes containing central black holes could result 
in the formation of a single super-massive spinning black hole — the preferred 
origin for powerful jets via the Blandford-Znajek mechanism (1977) — as well as 
the formation of the core of a luminous elliptical galaxy via heating associated 
with the orbital decay of the binary black hole system (Faber et al. 1997). 

Whatever the underlying physics, the evolutionary consequence of the nor- 
malisations of the luminosity functions of Fig. 2 and quasar lifetimes ~ 10^ yr 
is that every massive (^L*) galaxy was plausibly a quasar at high redshift (e.g. 
Richstone et al. 1998), and that ~ 1 in 100 of such galaxies spent time as a 



^ This 100:1 ratio of radio-quiet to radio-loud objects is somewhat higher than most other 
estimates in the literature (e.g. Goldschmidt et al. 1999). The main reasons for this are asso- 
ciated with the effects of Doppler boosting. First, both the optical and (high-frequency) radio 
emissions from an intrinsically radio-loud object can be boosted such that more common (but 
intrinsically lower Lboi) objects are preferentially promoted into an optically-selected sample of 
quasars if they have powerful relativistic jets. Second, low-power relativistic jets seem common 
amongst the radio-quiet population (e.g. Miller, Rawlings & Saunders 1993); when aligned 
with the line-of-sight, these jets produce sufficient radio emission that they can be designated 
as 'radio-loud' despite the absence of an intrinsically powerful jet. 
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powerful radio source. The present-day (co-moving) space density of the sites 
of past Eddington-tuned jet activity is given roughly by the density at 2; ~ 2.5 
times the ratio of the duration of 'the quasar epoch' (see Sec. 4) to the duration 
of a radio source outburst, or ~ 10~'^ X 4 X 10^/10*^ « 4 X 10~6 Mpc-3 which 
turns out to be precisely the space density of rich clusters (Dalton et al. 1994). 
Plausibly, therefore, every rich cluster of galaxies was once the site of extreme 
radio activity. 

4. Cosmic evolution at z > 2.5: the redshift cut-off 



311 ' ' ' ' ' I ' ' ' ' ' ' ' I 1.0 




Figure 3. (Left) Rough constraints on the high-rcdshift evolution 
of the luminosity density H = / i/Lj, [dp/d log j^q -^i^]d(logio Ly)di' of the 
Eddington-tuned radio population (which dominates the total radio lu- 
minosity density at high z). These have been calculated from model C 
of Willott et al. (2000c), amended at high-redshift to match the ~ 90 
per cent confidence intervals on the high- 2; decline in 11 inferred (indi- 
rectly) from constraints on the decline in the apparent space density 
of the flat-spectrum population (see Jarvis &: Rawlings 2000) — direct 
constraints (e.g. Jarvis et al. 2000) are even looser! The limits taken 
for the frequency integral were 10 MHz to 100 GHz in the rest-frame, 
and a = 1 was adopted. (Right) A rough calculation of the heat in- 
jected into the intracluster medium (ICM) by Eddington-tuned radio 
sources per co-moving Mpc~^ as a function of cosmic time. The solid 
line shows the prediction made on the basis of the lower bound to the 
high-z RLE in the left-hand figure, and the dotted line represents the 
upper bound. This calculation makes the following assumptions: (i) 
that the total jet power Q per Mpc^ can be obtained from the RLE 
using the transformation Q = 10^^'^ x L^§[' (Willott et al. 1999); and 
(ii) that 50 per cent of the jet power is dissipated as heat in the bow- 
shock of the radio source (Rawlings &; Saunders 1991; Willott et al. 
1999; Blundell & Rawlings 2000). 
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The results of Sec. 3 beg the question of constraints on the evolution of the 
radio source population at high redshift, i.e. at z > 2.5. I will focus my attention 
on the Eddington-tuned population since, as is obvious from Fig. 1, the lower- 
luminosity population is not detectable at high redshifts in bright radio-selected 
samples. My contention here is that, even though many radio-luminous objects 
have now been observed in the z > 2.5 Universe, constraints on their space 
density remain very loose. Since these points are discussed in detail elsewhere 
(Jarvis & Rawlings 2000; Jarvis et al. 2000; Willott et al. 2000c), I wiU simply 
summarise my view of the current situation. 

• Direct studies of the steep-spectrum population at z > 2.5 based on red- 
shift surveys of complete samples (Willott et al. 2000c) do not allow us 
to discriminate between models in which the T of radio sources declines 
at high redshift from those in which it stays constant. This is a subtly 
different conclusion to that reached by Dunlop & Peacock (1990), and the 
possible reasons for this difference are discussed by Willott et al. (2000c). 

• Through finding steep-spectrum radio galaxies at z ~ 4 — 5 (e.g. Rawlings 
et al. 1996), direct studies of the z > 2.5 population using 'filteredj^samples 
suggest roughly constant T from 2;~2.5to2;~5, although gradual 
declines or even inclines cannot be ruled out (Jarvis et al. 2000). 

• Indirect^vidence for a high-z cut-off in T from the study of luminous 
flat-spectrum radio quasars is more equivocal than received wisdom might 
suggest. Gradual declines in the space density of the most luminous flat- 
spectrum quasars, by a factor ~ 4 between z ~ 2.5 and z ~ 5 (Dunlop 
& Peacock 1990), are favoured by the data but the observable co-moving 
volume is too small to make deflnitive statements: both constant space 
density models and abrupt cut-offs, like those suggested by Shaver et al. 
(1996), are ruled out only at the w 2a level (Jarvis & Rawlings 2000). 

This state of ignorance concerning the high-z evolution of radio sources is 
illustrated in Fig. 3. It is worth noting that it is not clear that any other pub- 
lished measures of the high-z evolution of quasars are any better constrained, 
largely because of the difficulty of correcting for obscuration in optical/soft X- 
ray selected samples. In fact the abrupt cut-off found by optical quasar surveys 
(e.g. Warren, Hewett &; Osmer 1994) and the roughly constant space densities 
found by X-ray surveys (e.g. Miyaji et al. 2000) are together suggestive of ob- 
scuration: as a spectrum becomes increasingly redshifted, attenuation of optical 



A filtered sample is one in which optical spectroscopic follow-up is confined to a small fraction 
of the radio sources in a fiux-density limited sample, where the filters used by current surveys 
are a (steep) spectral index criterion and/or a (small) angular size criterion (e.g. Blundell et 
al. 1998). 

* Eddington-tuned jets pointed within a few degrees of the line-of-sight are Doppler-boosted 
to produce flat-spectrum radio quasars; these are thought to be the dominant fiat-spectrum 
population in high-frequency selected radio samples [although see Jarvis & Rawlings (2000) 
for a discussion of the importance of GHz-Peaked sources]. The apparent space density of 
these fiat-spectrum sources gives an indirect measure of the space density of the parent steep- 
spectrum population close to, but above, the break in the RLF (Jarvis & Rawlings 2000). 
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light becomes more effective due to tfie increased efficiency of dust obscuration 
at UV wavelengtfis, whereas the opposite is true for the effects of redshifting on 
the opacity at X-ray wavelengths (e.g. Wilman &: Fabian 1999). 

A large number of programmes are compromised by ignorance concerning 
the redshift cut-off, including basic questions like: (i) whether or not the onset 
of quasar activity re- ionized the Universe (e.g. Haiman & Loeb 1997); (ii) un- 
derstanding the nature of the sub-mm source counts (e.g. Almaini, Lawrence &i 
Boyle 1999); and (iii) using the density of high-z objects to constrain the power 
spectrum of density fluctuations on ~Mpc scales (e.g. Efstathiou & Rees 1988). 

I concentrate here, however, on an important physical parameter connected 
directly to the high-z evolution in radio luminosity density — the contribution 
of Eddington-tuned radio sources to the entropy budget of the Universe. The 
results of a rough calculation are shown in Fig. 3. The key point is that although 
the radio source phenomenon is 100— times rarer than the quasar phenomenon, 
the relativistic jets feeding radio sources have bulk powers which are of the same 
order as the Eddington luminosity iEdd of a ~ 10^ Mq black hole (Rawlings & 
Saunders 1991; Willott et al. 1999; Blundell & Rawlings 2000), and >100-times 
more powerful than the relativistic jets in luminous radio-quiet quasars (Miller, 
Rawlings & Saunders 1993). So, provided that the power in quasar winds is 
^0.01 i^Lcdd (Silk & Rees 1998) and provided one can neglect photoionisation 
heating (Valageas & Silk 2000) , it seems plausible that more power is available to 
heat the IGM as the result of radio source induced shocks, than is yielded by the 
more numerous radio-quiet quasar population. This heat, or entropy, input into 
the IGM will be concentrated in the clusters of galaxies which at one time hosted 
powerful radio sources. For the cut-off model in Fig. 3, the duration of the quasar 
epoch is 4 Gyr (~ 10^'^s) and a crude estimate of the heat input into each rich 
cluster by Eddington-tuned radio sources is ~ 10^^ x 10^^/4 x 10~^ 2.5 x 10^^ J. 
This is a reasonable fraction of the total thermal energy (~ 5 x 10^^ J) in a rich 
(~ 10^^ M0 in baryons) cluster. For the no-cutoff model of Fig. 3, the quasar 
epoch lasts longer, and the heat input approaches ~ 10 per cent of the total 
thermal energy in rich clusters. Understanding Fig. 3 is therefore likely to be 
central to any understanding of the entropy budget of the Universe. At earlier 
cosmic times the heat input by radio sources was probably sufficient to overcome 
the gravitational binding energy of gas in the relatively shallow potential wells 
of forming clusters; this would provide an attractive alternative to supernova- 
and quasar-driven mechanisms invoked to explain the 'excess' entropy (over and 
above that due to shock heating during cluster collapse) identified by studies of 
the properties of rich clusters (e.g. Kaiser 1991). In essence I am suggesting a 
feedback mechanism similar to the quasar-driven wind model invoked for galaxy 
formation by Silk &: Rees (1998) but applicable to more massive systems, i.e. 
rich clusters, and to more powerful outflows, i.e. radio source driven bow shocks. 

5. The sub-mm view of the cosmic evolution of radio sources 

The strong cosmic evolution of the Eddington-tuned quasar population discussed 
in Sec. 3 is most likely the simple and direct consequence of the standard hi- 
erarchical model for galaxy formation in which major mergers of galaxy-sized 
sub-clumps are more prevalent in the younger Universe. A prediction of this 
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Figure 4. (Left) L151 versus z plane for 47 radio galaxies [one at 
logio(-^^i5i) = 26, z = 1.43 is off the plot], and one radio quasar 3C318, 
observed at 350 GHz with SCUBA (triangles - 3CRR objects, circles 
- 6CE objects, squares - objects from other surveys); filled symbols 
represent detected objects, and open symbols those with upper limits. 
The sizes of the symbols are scaled from the 2a 350 GHz flux density 
5350 Hmit (^350 ~ 2 mJy) by the factor [1 + ln(S'35o/2)]. (Right) Illus- 
tration of the fraction of global star-formation associated directly with 
the quasar population. The star-formation associated with radio-loud 
objects is shown by the two lower lines adopting either the lower bound 
to the high-z RLF (solid line, the cut-off model of Fig. 3) or the upper 
bound (dotted line, no cut-off model) . This has been calculated by as- 
signing a star formation rate (SFR) of 10(1 -|- z)^ Mq yi"^ (saturating 
at 2; = 4) to each Eddington-tuned radio source, without any correla- 
tions between L151 and SFR. To account for the radio-quiet population 
the upper lines scale these functions by a factor of 100 (see Sec. 3) . The 
symbols represent the following: star, local star formation rate from 
the Ha survey of Gallego et al. (1995); circle, SCUBA HDF survey 
(Hughes et al. 1998), given by Hughes et al. as a lower limit but likely 
to be within a factor of a few of the correct value given the brightness 
of the sub-mm background; square, global SFR corresponding to one 
of the SCUBA HDF sources lying at 4 < 2 < 6. 



hypothesis is that the Eddington-tuned radio population is associated with sites 
of intense star formation since, at high z, a large fraction of the baryons in the 
sub-clumps are likely to be in the form of gas rather than stars (e.g. Storrie- 
Lombardi, McMahon & Irwin 1996). The observed properties of high- 2; radio 
sources should show signs of this. Although UV and optical observations have 
now provided at least one compelling example of a huge starburst associated 
with a high- 2; radio galaxy (Dey et al. 1997), it has proved impossible to prove 
the ubiquity of the starburst phenomenon from UV/optical/ncar-IR, or indeed 
radio, observations because any starburst signature can be so easily swamped 
by emission intimately associated with the active nucleus, e.g. in the case of UV 
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observations, scattered light from the quasar nucleus. A far cleaner diagnostic 
of large rates of star formation in radio galaxies is a huge rest-frame far-infrared 
luminosity^, requiring a survey of high- 2: radio galaxies in the sub-mm wave- 
band. The advent of the SCUBA camera on the JCMT has recently made such 
a survey possible. Full details of this work are published elsewhere (Archibald 
et al. 2000), so I concentrate here on the major results and their implications. 

Fig. 4 shows how the chances of making a sub-mm detection of a radio 
galaxy varies across the L151-Z plane. Considering z>3 and log;^o(-^i5i)~28.5, 
we see that all the radio galaxies are detected by SCUBA, whereas at lower 
redshifts and log^Q(Li5i) the fraction of detections falls abruptly to ^20 per 
cent. Archibald et al. make a detailed study of the underlying correlations 
between rest-frame far-IR luminosity LfiRj -^^isi and z. They conclude that the 
major effect at work is an ~ (1 + z)^^ rise in -LpiR for radio galaxies out to z>4. 
The most straightforward interpretation of this result is that the characteristic 
star formation rate (SFR) of luminous radio galaxies rises from £10 Mq yr~^ 
at 2; = to >1000 Mq yr"-*^ at z ~ 4, an interpretation which seems to fit in 
naturally with the arguments made in Sec. 3. If the strong cosmic evolution 
in T for the Eddington-tuned radio population is indeed driven by systematic 
changes in the trigger rate with z then the most viable trigger mechanisms, e.g. 
the formation of a rapidly spinning super-massive black hole in a major merger 
event, are likely to be far more prevalent at high z; provided, as seems likely, the 
progenitor galaxies are gas rich, intense associated starbursts seem inevitable. 

In Fig. 4 I take these arguments to their logical conclusion by deriving rough 
estimates of the star formation in the Universe associated directly with AGN. It 
seems that, as noted by Almaini et al. (1999), a large fraction (^10 per cent) of 
the star formation in the high-z Universe is probably occurring directly alongside 
quasar activity in the cores of forming galaxies (see also Boyle & Terlovich 1998). 
Indeed, if we hypothesise that galaxies form in a series of short (less than or of 
the order of a quasar lifetime^), intense (JiJlOO Mq yr~^) starbursts across the 
entire quasar epoch (a process which perhaps culminates in the quasar phase) 
then the important star formation at high redshift goes hand-in-hand with the 
formation of black holes, and the epoch of major quasar activity. Of course if 
Compton-thick objects turn out to be the dominant AGN population at high-z, 
then the link between AGN and star formation activity may prove closer still. 



^ Photometric measurements on the Rayleigh- Jeans side of a thermal spectrum can be used 
to confirm that the emission is due to dust (which for radio galaxies it almost always is) and 
provide good estimates of dust mass which, in well-studied objects, correlates closely with gas 
mass (e.g. Hughes, Dunlop & Rawlings 1997). Since it is often difficult to prove beyond doubt 
that the dust is heated by stars rather than an active nucleus (e.g. Willott, Rawlings & Jarvis 
2000e), the correlation between dust/gas mass and star formation rate is more uncertain. 

^ One of many caveats to my interpretation of Fig. 4 concerns the question of the relative 
timescales of AGN and star- formation activity. Blundell & Rawlings (1999) have pointed out 
that selection effects ensure that the known z > 3 radio galaxies, e.g. those in the SCUBA 
survey of Archibald et al. (2000), are viewed within 10^ yr of the jet-triggering event. If the 
characteristic timescale of the starburst activity is shorter than the ~ 10* yr timescale of AGN 
activity, then the 'youth-redshift degeneracy' (Blundell & Rawlings 1999) may be an important 
contributor to the observed 'cosmic evolution' in the SFR of powerful radio galaxies (Fig. 4). 
It is interesting in this regard that one of the youngest 3CRR radio sources at 2 ~ 1.5, 3C318, 
is a SCUBA-detected quasar (Willott et al. 2000e). 
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The results of this section concerning the history of star formation, and of 
Sec. 4 concerning the global entropy budget, should provide food for thought for 
the growing band of researchers using semi-analytic models to explain galaxy 
formation and evolution (e.g. Somerville &; Primack 1999). At present these 
models take no account of feedback between quasars, their jets and their host 
galaxies. For observers of high-2 AGN it is obviously vital to tie down huge and 
basic uncertainties like the evolutionary behaviour of quasars at z > 2.5. 

6. Future prospects: bigger redshift surveys and better modelling 

In previous sections we have seen how cosmological studies using low-frequency 
selected radio sources are compromised by the limitations of the current gener- 
ation of redshift surveys. With these considerations partly in mind Gary Hill 
(TEXas) and I (Oxford) have begun a programme of TexOx redshift surveys. 
One aspect of this programme — the TexOx-1000 (or TOOT) survey — is to 
measure ~ 1000 redshifts for a complete sample flux-limited at >S'i5i = 100 mJy. 
The optical/radio/near-IR imaging part of this survey is now almost complete 
and spectroscopic follow-up began in earnest in January 2000. This programme 
aims to exploit the complementary capabilities of UK facilities like the ISIS 
spectrometer on the WHT - the instrumental back-bone of the 6CE and 7CRS 
programmes - and the LRS spectrometer on the new 9.2-m Hobby-Ebcrly Tele- 
scope. We are already about a quarter-way through the required spectroscopic 
follow-up, and we anticipate completion of the TOOT survey in mid 2002. As 
illustrated by Fig. 1, the TOOT survey will be the first to directly measure the 
space density of logiQ(Li5i)>26.5 objects, and hence the luminosity density of 
the steep-spectrum population, at z ~ 2.5. 

Still larger redshift surveys may be needed to settle the question of evolution 
at z ^ 2.5. A major difficulty of making such surveys is the lack of any multiplex 
capability at relatively bright flux-density limits. At the 5i5i ~ 100 mJy limit of 
the TOOT survey, the average separation of sources on the sky is ~ 20 arcmin so 
there is still no effective alternative to object-by-object spectroscopic follow-up, 
particularly as the largest field-of-view multi-object systems tend to compromise 
one or more of sensitivity and spectral coverage. In the medium-term future, 
progress can probably best be made by refining photomctric-rcdshift techniques 
for the radio source population, although the requirement for wide-field capabil- 
ity at near-IR wavelengths means that substantial progress will probably require 
new initiatives like the planned surveys with VISTA (see Jarvis et al. 2000). In 
the long-term future, the hope is that redshifts will come for free as facilities like 
the Square Kilometre Array push radio sensitivities to the level where surveys 
can routinely detect HI and CO features superimposed on the radio continuum. 

As we await the results of bigger redshifts surveys, a high priority is to 
build better models for the cosmic evolution of the radio source population. 
This splits into two categories of problem. First, we need to develop robust 
models for the physical evolution of radio sources: i.e. once a source of a given 
jet power has been triggered, how does its radio properties (e.g. luminosity, 
size and spectral index) develop as a function of time, and as a function of the 
radio source environment? Second, we need to understand systematic changes 
in radio source environments with z, both because of their influence on radio 
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properties, but more importantly so that we can understand how jet trigger rates 
reflect the build up of structure in the Universe. In the first area, the hiatus 
following the pioneering work of Scheuer (1974) is now well and truly over (e.g. 
Kaiser, Dennett-Thorpe & Alexander 1998; Blundell et al. 1999). In the second 
area, semi-analytic models could provide the way forward, so these must now 
be developed to incorporate the feedback between AGN, their jets and their 
environment. 
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